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Выделение скрытых областей 
слабоконтрастных изображений 
методом ортогональной декомпозиции 


Рассмотрены информационные возможности метода анализа слабоконтрастных изображений. Суть 
метода основана на разделении анализируемых компонент многопараметрового (мультиспектрального) 
изображения на ортогональные составляющие, при этом в составляющих высокого порядка выделяются 
области, первоначально неразличимые на исходных изображениях. Представлены результаты численного 
моделирования и реальной проверки работоспособности метода. 


Введение 


Задача выделения визуально неразличимых областей слабоконтрастных изобра- 
жений (как обычных, так и многопараметровых) имеет большое прикладное значение 
в таких областях, как медицинская диагностика (идентификация радиологических, 
ультразвуковых и оптических изображений), дистанционное зондирование Земной по- 
верхности и анализ изображений геофизических полей. Наличие подобной проблемы 
обусловлено недостаточной чувствительностью человеческого глаза, воспринимаю- 
щего, в соответствии с законом Вебера, лишь 2% перепадов яркости, причем это значе- 
ние существенно зависит от окружающего фона, причем сложность задачи существенно 
возрастает с увеличением числа диапазонов (параметров) в мультиспектральных изоб- 
ражениях. 

Целью работы является демонстрация информационных возможностей метода 
ортогональной декомпозиции применительно к решению задач выделения визуально 
неразличимых участков слабоконтрастных изображений. 


Ортогонализирующие преобразования 


Для определенности первоначально рассмотрим случай, связанный с анализом 
многопараметрового (например, цветного) изображения { Г(х, у, п) }, где п- число па- 


раметров (диапазонов), соответствующих каждому пикселю (х, у). Если представить 
анализируемый ансамбль изображений {1} в виде прямоугольной матрицы А размер- 
ностью (М хМ№)хи, где Ми М№- число пикселей вдоль осей «х» и «у» соответствен- 
но, то можно произвести сингулярное разложение матрицы А на основе выражения 
<=и\ЩУу`, (1) 
где 0 = [61 те, я — матрица сингулярных векторов размерностью ((М * №)хп) , фор- 
мирующих ортонормированный базис пространства, натянутого на столбцы матрицы С, 
формирующая базис «собственных изображений» (СИ) [1]; И’ - диагональная мат- 
рица, элементы которой и’ (1=1.....п) являются сингулярными числами матрицы С, 
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причем у, 2 м, >...> у, >20. Они являются собственными значениями матрицы С”С, 


а Г представляет собой ортогональную матрицу размерностью (пхи). На основе ана- 
^ 


лиза спектра нормированных сингулярных чисел и; базируются основные методы 
компрессии (сокращения избыточности) ансамбля анализируемых изображений {1}, 
основанные на том обстоятельстве, что если 


ии. С) 
1=1 


где О - критическое значение, равное, например, 95% или 99%, А < п. Физический 
смысл выражения (2) базируется на том обстоятельстве [2], что «собственные изоб- 
ражения», синтезированные из столбцов матрицы И (т.е. Ц`,....И‚), обеспечивают 


р% всей информации, заложенной в исходном ансамбле {П. 
На рис. 1 представлена численная модель трехдиапазонного изображения. 


МИР 


Рисунок 1 — Численная модель трехдиапазонного изображения 


Характерной особенностью рассматриваемого примера является то обстоятель- 
ство, что на третьем изображении расположено семь слабоконтрастных участков, из 
которых визуально идентифицируются только три, а четыре являются невидимыми. 
На рис. 2 представлены «собственные изображения» анализируемого ансамбля, спектр 


нормированных сингулярных чисел которого равен: и=85,3405%; и›=14,6582%; 


из =0,0012%. Из рассмотрения спектра сингулярных чисел следует, что первое «СИ» 
содержит более 99% информации, содержащейся в трех исходных изображениях, т.е. 
в соответствии с рекомендуемым теоретическим подходом (2) можно ограничиться 
только анализом (.. 


1 2 3 


Рисунок 2 — «Собственные изображения» ансамбля, представленного на рис. 1 
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Однако из рассмотрения рис. 2 видно, что визуальная идентификация всех семи сла- 
боконтрастных участков возможна лишь на третьем СИ, содержащем лишь 0,0002% 
исходной информации. Из рассмотренного примера вытекают два практических важ- 
ных вывода. 

1. Для обнаружения и идентификации слабоконтрастных и визуально неиденти- 
фицируемых участков необходимо анализировать весь ансамбль синтезированных СИ. 

2. Для информационной значимости ансамбля СИ необходимо, чтобы исходные 
изображения были высококачественными, поскольку в противном случае большая часть 
СИ будет соответствовать шумовым компонентам [3], не поддающимся анализу. 

Метод сингулярного разложения (1) обеспечивает алгебраическую ортогонализа- 
цию. Возможен и другой, статистический, подход, базирующийся на синтезе ансамбля 
«независимых изображений» — метод независимых компонент (МНК) [4]. 

МНК - это метод нахождения скрытых факторов или компонент в многомерных 
данных, особенность которого в том, что он позволяет выделять статистически неза- 
висимые составляющие. В принципе в классическом факторном анализе [5] также 
требуется, чтобы выделенные факторы были статистически независимы. Однако при 
этом предполагается, что анализируемые данные имеют гауссово распределение, ко- 
торое редко встречается на практике. 

Основной принцип нахождения независимых компонент базируется на максими- 
зации негауссовости анализируемого статистического распределения. Из центральной 
предельной теоремы следует, что сумма независимых случайных переменных стре- 
мится к нормальному распределению. Тогда сумма даже двух независимых случайных 
переменных обычно будет иметь распределение более близкое к нормальному, чем 
распределение любой из двух исходных случайных последовательностей. Задачу обыч- 
но называют задачей слепого разделения источников, подчеркивая, что единственной 
информацией для восстановления входа является реализация вектора наблюденийх. 

В МНК предполагается, что анализируемые случайные переменные х,....,х, могут 


быть представлены в виде линейной комбинации случайных переменных 3+....,5, как 
х,=а15, +а,.5. +...а5, и =12,..тЕ<п, (3) 
где а, - некоторые неизвестные нам действительные коэффициенты; 5, — рассмат- 


риваются как статистически взаимно независимые скрытые компоненты. В векторно- 
матричной форме записи уравнение (3) представимо в виде х = 45. 

Основные допущения МНК базируются на трех основных предположениях. 

1. Независимые компоненты предполагаются статистически независимыми, Т.е. 
совместная плотность распределения вероятностей р(5,5.,...,5,) может быть факто- 
ризована в виде 

Р(51,55,....5%) = р($.)р($>)...р($ь). 

2. Независимые компоненты должны иметь негауссово распределение и прибли- 
зительно идентичные распределения. Поэтому если искомые компоненты имеют га- 
уссово распределение, то реализация МНК невозможна. 

3. Неизвестная матрица А может быть инвертирована, т.е. возможно обращение 
уравнения (3) 

В и. 

Первое из предположений является гораздо более мощным по сравнению с тре- 
бованием некоррелированности (как в методе главных компонент). Некоррелирован- 
ность (в случае произвольного распределения плотностей вероятностей) не означает 
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статистической независимости, в то время как если, например, две переменные $1 и 5, 
статистически независимы, то любые их нелинейные преобразования типа 2(5,) и й(5,) 


будут одновременно и некоррелированными. 
Для данного вектора х задача сводится к поиску разделяющей матрицы В ‚ такой, 
что исходный вектор 5 может быть получен из выходного вектора у (рис. 3). Это зна- 


чит, что решение можно представить в виде у= Вх = ВА5. Если обозначить одну из 


строк матрицы В как Ь, то линейная комбинация Б'х будет эквивалентна одной из 
независимых компонент. 


Вектор Вектор 
источника наблюдения 
5 х 


Выходной 
вектор У 


Смеситель 


А 


Разделитель 


В 


Рисунок 3 — Схема метода независимых компонент 
(«слепого разделения источников») 


Проблема состоит в том, как воспользоваться выводами центральной предель- 
ной теоремы для нахождения Ь? На практике мы не можем точно определить вектор 
Ь, так как нам ничего не известно о матрице А. Все упирается в необходимость на- 
хождения аппроксимационного решения, заключающегося в выборе такого вектора 
Ь ‚ который бы максимизировал негауссовость распределения Б'х ‚ дающего одну из 
независимых компонент. 

Обычно в этой задаче предполагается, что входной сигнал имеет нулевое сред- 
нее значение, что гарантирует нулевые средние для хи выходного вектора. 

Классической мерой негауссовости является кумулянт 4-го порядка, или эксцесс 


у, =Е“}- ЗЕ”, (4) 
где Ё 4х" — п-й математический момент. 


Однако он является весьма чувствительным к влиянию шумовых факторов. По 
этой причине в экспериментальных исследованиях, в качестве меры негауссовости, 
можно использовать значение энтропии. Энтропия случайного вектора 5 с плотностью 
распределения р(5) определяется как 


Н($) =-| рю рб). (5) 


Значение энтропии для МНК обусловлено фундаментальным выводом теории 
информации - гауссовы переменные имеют наибольшую энтропию (т.е. наибольшую 
непредсказуемость) среди всех случайных переменных с эквивалентными дисперсия- 
ми. Практически именно это обстоятельство и означает, что значение энтропии мо- 
жет быть использовано в качестве меры негауссовости. 

Определение матрицы В достигается путем максимизации значения отрицатель- 
ной энтропии Н(5) в предположении, что 5 = / (5), где Г рассматривается как ап- 


проксимирующая функция распределения независимых компонент. 
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Таким образом, структура использованного алгоритма МНК включает в себя сле- 
дующие этапы. 

1. Центрирование исходных данных для обеспечения нулевых средних. 

2. Проведение предварительного декоррелирующего преобразования с определе- 
нием значимого числа главных компонент, что позволяет минимизировать влияние 
шумовых факторов для последующего анализа (декоррелирующее преобразование вхо- 
дит как составная часть в МНК). 

3. Инициализация начального значения матрицы Ву случайными значениями с 


единичной дисперсией. 

4. Максимизация значения неотрицательной энтропии методом сопряженного гра- 
диента — в результате итерационного процесса находится оптимальная, с точки зре- 
ния критериев МНК, матрица В. 

5. Нахождение значения независимых компонент на основании уравнения (3), чис- 
ло которых совпадает с числом использованных главных компонент. 

На рис. 4 представлены три «независимых изображения», синтезированных из ис- 
ходного ансамбля на рис. 1. 


Рисунок 4 — «Независимые изображения» ансамбля на рис. 1, синтезированные 
на основе метода независимых компонент 
Из рассмотрения рис. 4 видно, что в отличие от метода сингулярного разложе- 
ния (рис. 2) «аномальные» участки выделились на всех трех синтезированных изоб- 
ражениях. На рис. 5 представлены результаты ортогональной декомпозиции цвет- 
ного ВОВ дерматологического изображения (рак кожи). 


иг 5 


Рисунок 5 — Дерматологическое цветное изображение: верхний ряд — 
КОВ компоненты оригинала; нижний ряд — синтезированные методом 
независимых компонент новые изображения 
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Здесь, как и в случае цифровой модели при использовании метода сингулярного 
разложения, «скрытая» информация выделилась на третьем синтезированном изобра- 
жении. Его диагностическая значимость остается под вопросом, поскольку результат 
явился неожиданным для самих врачей. 

Тот факт, что для «тонкого» визуального анализа необходимо использовать весь 
ансамбль синтезированных «собственных» или «независимых» изображений, откры- 
вает возможность повышения чувствительности визуального анализа и обычных (од- 
нопараметровых) слабоконтрастных изображений, базируясь на идее их автоморфного 
отображения [6], [7]. 

На рис. 6 а) представлено изображение оптического излучения из конца оптово- 
локонного кабеля. 


а) 6) 
Рисунок 6 — Изображение оптического излучения из конца оптоволоконного кабеля: 


а) оригинал; 6) второе «собственное изображение» автоморфного отображения 
оригинала 


Из рассмотрения рис. 6 6) видно, что на втором «собственном изображении» вы- 
делились особенности топологии рис. 6 а), а также скрытая деталь, расположенная в 
нижнем правом углу изображения. 

На рис. 7 а) представлено инфракрасное изображение работающей интегральной 
микросхемы. 


а) 
Рисунок 7 — Инфракрасное изображение работающей микросхемы: 
а) оригинал; 6) третье «независимое изображение» 
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Из рассмотрения рис. 7 6) видно, что, как и в случае оптического изображения 
на рис. 6, особенности анализируемого изображения выделились на третьем «независи- 
мом изображении» автоморфного отображения рис. 7 а). В данном случае применение 
как МНК, так и сингулярного разложения показало идентичные результаты. 


Заключение 


На основании проведенных исследований можно сделать следующие выводы. 

1. Применение метода ортогональной декомпозиции как многомерного ансамбля, 
так и обычного (однопараметрового) слабоконтрастного изображения позволяет вы- 
делять визуально неразличимые скрытые участки. 

2. Метод сингулярного разложения и метод независимых компонент можно счи- 
тать не конкурирующими, а дополняющими друг друга, поскольку метод сингуляр- 
ного разложения является намного более простым в вычислительном отношении по 
сравнению с методом независимых компонент. 

3. Применение метода ортогональной декомпозиции предъявляет повышенные 
требования к качеству (помехоустойчивость и точность регистрации) исходных изоб- 
ражений. 


Литература 


1. Форсайт Д. Компьютерное зрение: современный подход / Д. Форсайт, Ж. Понс ; [пер. с англ. 
А.В. Назаренко, И. Ю. Дорошенко]. — М. ; СПб. ; К : Вильямс, 2004. — 926 с. 

2. Фукунага К. Введение в статистическую теорию распознавания образов / Фукунага К. ; [пер. И.Ш. Тор- 
говицкого] ; под. ред. Ф.Ф. Дорофеюка. - М. : Наука, 1979. — 367 с. 

3. Ахметшина Л.Г. Адаптивная фильтрация шумов в сигналах и изображениях: метод селективного 
сингулярного разложения автоморфного отображения / Л.Г. Ахметшина // Искусственный интел- 
лект. — 2005. — № 3. - С. 328-335. 

4. Нууагтеп А. шдерепдеп! Сотропепе Апа1уз1$ / Нууагпеп А., Кафипеп /., Ола Е. — М.У. : лова \Пеу 

& 501$ Шшс., 2001.-354 р. 

Иберла К. Факторный анализ / Иберла К. ; [пер. В.М. Ивановой]. - М. : Статистка, 1980. - 398 с. 

6. Ахметшина Л.Г. Повышение разрешающей способности изображений геофизических полей на ос- 
нове метода многомерной ортогональной адаптивной кластеризации / Л.Г. Ахметшина // Науковий 
Висник Нацюнального грничого унверситету. - 2003. — № 10. - С. 35-38. 

7. Ахметшина Л.Г. Повышение чувствительности анализа рентгеновских снимков методом автоморф- 
ного отображения в базисе собственных изображений / А.М. Ахметшин, Л.Г. Ахметшина // Кли- 
ническая информатика и телемедицина. - 2008. - Т. 4, № 5. - С. 30-36. 


— 


О.М. Ахметшин, К.О. Ахметшин 

Вид1лення прихованих Длянок низькоконтрастних зображень 

методом ортогонально? декомпозицИ 

Розглянуто 1нформащйн! можливост! нового методу анал1зу низькоконтрастних зображень. Суть методу 
полягае у роздленн! аналзованих компонент багатопараметрового зображення на ортогональн! складовт, 
при цьому на складових @лышп високого порядку видмляються длянки, невидим! на первсних зображеннях. 
Представлен! результати реально! перев1рки працездатност! методу. 


А.М. Акптебйуп, К.А. АКитебйуп 

Н1адеп Ротат$ Ехгасйоп оЁ Го\ Сопыга$ Ппасез Бу Меап$ о? Ог®огопа! Эесотро$#оп Мето4 
шЕЮюгтайоп роз Ш@ез оРа пе\и ше#о4 10% сопёга{ птасез апа|уз1$ аге соп$1еге4. ТВе тат 14еа оЁ фе 
тефоа 15 Юоипде4 оп зерагайоп о# апайуте ппазе оп огфогопа] сотропеп®. Рог зо[айоп оф {515 ргоМет 
ш Фе агаее ууеге изе4 Ко арргоасПез: Фе Ягз( опе 15 Базе оп изшя фе зтешаг уаше десотроз!Шоп (1 1$ 
ап а[еебга1с тефо4) ап Фе зесоп4 опе 15 Базе4 оп ш4ереп4еп{ сотропеп тео4. ТВе 1а5ё ргоу14ез Фе 
эюсвазйс огрогопа|тайоп. Н144еп дотайл$ о# апа[у7е4 ппазе аге зерагайе4 оп огфогопа[уте4 итазе оЁ Шов 
отаег. ТБе гези5 оё теа] {езйп® оЁ фе тефо4 роз Шйез аге ргезещед. ТЬе опез Ва зво\уп фа Бой тефо4$ 
аге сотр!етепЕ еасп о'ег. 
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